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We study a class of models for tri-bimaximal neutrino mixing in 50(10) grand unified SUSY 
framework. Neutrino masses arise from both type-1 and type-II seesaw mechanisms. We use dimen- 
sion five operators in order to not spoil tri-bimaximal mixing by means of type-I contribution in the 
neutrino sector. We show that it is possible to fit all fermion masses and mixings including also the 
recent T2K result as deviation from the tri-bimaximal. 
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I. INTRODUCTION 



Neutrino mixing leads to large atmospheric angle (maximal) , large solar angle (trimaximal) and small reactor angle. 
In particular recently T2K collaboration has given indication of non zero reactor angle 1}. After such a result, the 
. global fits of neutrino parameters give non zero reactor angle at 3cr: 



sin^ = 0.013™ sin2 0i3 = O.O25t°f2^ i. (1) 



This interesting result seems in contradiction with tri-bimaximal (TBM) mixing ansatz [J| that predicts zero reactor 
angle. However TBM solar and atmospheric mixing angles can be used as first approximation. Deviation from zero 
reactor angle can arises in grand unified theory (GUT) like SU{5) and 5*0(10) from the charged sector, see for 
instance [5] . While neutrino mass matrix is diagonalized from TBM unitary matrix, charged leptons are not diagonal 
giving deviation to the TBM. In this paper we consider such a possibility in the framework of a supersymmetric 
(SUSY) S'O(IO) model. In this scenario charged leptons and CKM mixings are strongly related, we therefore consider 
, the TBM as a good starting point to be corrected in general by small (CKM-like) deviations. In Ref. 6] has been 
CO ■ shown that in a renormalizable S'O(IO) model this is not possible in case of type-I seesaw. Such a difficulty arises from 
the fact that up quark and Dirac neutrino Yukawa couplings are strongly related in renormalizable 30(10) models. 
Some interesting attempts to obtain TBM with a flavour symmetry are developed in Ref.Q and Ref.Q, assuming 
T— I [ type-II seesaw to be dominant ^. In Ref.[l^, in the context of SUSY renormalizable iS'O(lO) with type-II seesaw 
T-H . dominance, a fit of all the fermion masses and mixing has been done (see also [isj). The superpotential considered is 
^ ' of the form 
• ^ . 

^ ■ w = /il6 16 10 -I- /16 16 126 + h'l6 16 120, (2) 

5t where /i is a symmetric matrix, h' is antisymmetric and / has the TBM structure, namely 

/2 /i fi 

/ = I A /2 + /O /2 - /O I , (3) 

/i /2 - fa f2 + fo 

where m^i = A ^ /i, "^1^2 = ./2 + 2/i and 771^,3 = /2 — /i + 2/o. It is well know that a mass matrix with the above 



structure is diagonalized by TBM mixing matrix, see for instance |14| . No assumptions have been made taking / to 
be TBM because we can always go to this basis by rotating the 16 of fermions [12!|. The matrices h, h' are assumed 
to be hermitian that can correspond to an underlying parity [l5l |. 



'Electronic address: ''blankenburg@fis .uniroma3.it | 
^Electronic address: morisiOific.uv.esl 

^ For an incomplete list of papers with TBM in GUT see Q|, for 50(10) models with discrete flavour symmetry and no TBM mixing see 
[lOl | and for a collection of general 50(10) models see at the References in Ref. [TH 
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Another possibility to reproduce TBM mixing in the framework of 5*0(10) GUT models is to use non-renormalizable 
operators containing a scalar field transforming as a 45// of 50(10) [y, [la, [13 • This field allows to distinguish up 
quarks from neutrinos permitting TBM mixing also in the case where neutrino masses arise from type-I seesaw 
mechanism. In particular in Ref. 6] for this purpose the dimension five operator 16 16 120/f 45// has been used. This 
operator yields a contribution to the up-quark mass matrix and not to the Dirac neutrino one allowing to distinguish 
the up-quark from Dirac neutrino sectors. In this way it is possible to obtain both Dirac and Majorana neutrino 
masses TBM and hierarchical structure in charged fermions sector. 

A full fit of quark and lepton masses and mixing in models with TBM mixing from type-I seesaw in 5*0(10) is 
still missing. In this paper we consider such a problem. We link the idea of distinguish up-quark and Dirac neutrino 
by means of 16 16 120// 45// operator with the result of Ref. 12| where has been shown that from the superpotential 
([2]) it is possible to fit all the data having TBM mixing in the neutrino sector. In this paper we will translate the 
superpotential © in the language of dimension five operators. 



In the next section we will review some of the 5*0(10) dimension five operator that will be useful to construct 
an 5*0(10) model giving TBM mixing with type-I as well as type-II seesaw following the indication given in the 
superpotential © of Ref. pjj]. In section lllll we give some examples of models and the corresponding fits, in sectionlTVl 
we discuss the possibility to obtain a renormalizable model, then in sectionlVl we give our conclusions. 



II. DIMENSION FIVE EFFECTIVE OPERATORS 



In this section we report the main ingredients that will be useful to construct our model in the next section. It 
contains some of the result of table VIII of Ref.[^ that we report in appendix A for the useful of the reader. 

As discussed in the introduction, one possibility to reproduce TBM mixing in 50(10) in the case of type-I seesaw 
is by means of the dimension five operator 16 16 120// 45// that allows to distinguish between up-quark and Dirac 
neutrino sectors. In this section we remark the feature of some dimension five operators that we will use in the next 
section. 

In general an 50(10) dimension five operator can be written as 16 16 (/)q (?iif, where (t>a,b sue scalar fields (j)a.b = 
1//, 16//, 16//, 45//,... and so on. For simplicity we assume that 50(10) is broken through SU{5) and we describe 
the contribution of the dimension five operators to the fermion mass matrices in the SU{5) language. When one of 
the components of ipa and (pi, take vev a, and bi respectively (where i is the 5*C/(5) index of the component), one 
generates contributions to the quark and lepton masses. Note that and bi can be possibly equal if (j)a — (j^b- 



The dimension five operators that will be used are: 

• 1616T6//T6// 

This operator can be obtained, for example, integrating out a 50(10) singlet 1^ or a 45-plct 45^ of heavy 
messenger fermions: 

/(16T6//)i(16T6//)i, /(16T6//)45(16T6//)45. (4) 

^From the table in appendix A it is possible to see that the first operator contributes to Y^, M/j and Ml while 
the second one contributes to Y^, M^, Ml and Yu- 

We assume / to have the TBM form(l3|), in this way the light neutrino mass is TBM with type-I or type-II 
seesaw, see also eq. ([M)) below. Note that / TBM is general because we can always go in this basis by rotating 
the 16 

• 1616120//45// 

This operator can be obtained by integrating out a couple 16^ — 16^ 

/i(16120//)i6(16 45//)t6. (5) 

It can yield a contribution to the up-quark mass matrix (and to the down-quark and charged lepton mass 
matrices) and not to the Dirac neutrino one, allowing to distinguish the up-quark from Dirac neutrino sectors. 

This can be described naively in the SU{5) language as follows. The up-type Higgs doublet in a 45// of SU{5) 
(contained into the 120//) couples antisymmetrically to the two matter multiplets 10 that give the up-quark mass 
Mu while it does not contribute to the Dirac neutrinos, when contracted as in eq. ([5]). This is a well know feature 
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of 5*0(10) where the operator 16 16 120^ contributes to cx {4:5su(5)) and cx {5su{5)) with antisymmetric 
Yukawa. But with the dimension five operator the resulting mass matrix is not antisymmetric and so can give 
the hierarchical structure of the quark sector. In fact with the insertion of a 24^f scalar multiplet (contained in 
the 45// of 5*0(10) ) that takes vev in the hypercharge direction, results that Af„ is not antisymmentric (but a 
generic matrix) since the Clebsch-Gordan coefficients for the isospin doublet Q and the isosinglet u'^ are different 
with a relative factor (—4) that arise from the hypercharge. 

We can describe the property of the 16 16 120/f 45// operator in more detail as follow. The 45// can take vev in 
its singlet lsu{5} component called X-direction^ or along the adjoint 24su{5) component, that is the hypercharge 
F-direction (see for instance (3). We indicate their vevs as 

bl = {^SU(5)), ^24 = (245[/(5)). (6) 

Equivalently, the 45// can take vev along the isospin direction or the B — L direction and their corresponding 
vev are denoted as 63 and 615 respectively and are given by 

&i = i(fo3 + 36i5) , 624 = J(-63 + 2fei5) . (7) 

The SU{5) components of the 120// of 50(10) that contain SU{2) doublet (giving rise to the Dirac masses 
terms for the fcrmions) are the 4,5su{5)i 455/^(5), 5su(5) ^^^d 5su{5) representations. We denote their vevs as 

05 = (5sc/(5)), 05 = (5sc/(5)), 045 = (455//(5)), 045 = (45sc/(5))- (8) 
iFvom the table in appendix A we have for instance that 

Mu = /l 045 (61 -4^24) - /l^ 045(61 +&24), (9) 

A'U = 5/10561 - /i^ 05 (-3&1 - 3&24), (10) 
Md = h{a^+ajg){~3bi+2b24)-h'^{a^+aj^)ibi+b2i), (11) 

MJ ^ h{a^^'Sa-^){-3bi-3b24)-h^ia-g^3ajg){bi+6b24). (12) 

Then if 05 = this operator contributes to Y^, Yd, Y^ and not to Y^. So h can be in a hierarchical form, with 
(3, 3) element dominant, as required by charged fermion phenomenology, without changing the TBM result. 
Because of the 45//, the mass matrix that results from such an operator is not antisymmetric but general. Note 
that with a 45// in B-L direction the resulting mass matrix is symmetric, in fact using eq. ^ we have 

Mu = /i 045 (63 - 615) - ^^045(615), (13) 

= /ia5(&3 + 36i5) + /i^a5(36i5), (14) 

Md = h{a-g+ ajg){-b3 - bis) - h'^{a-g+ ajg)bi5, (15) 

Mj = hia^~3aj^)i~3bi5)~h^ia^~3ajg)i3bi5~b3). (16) 

and stetting 63 = then all the mass matrices are symmetric. 

• 161610//45// 

This operator can be obtained by contracting with a couple 16^ — 16^: 

/i'(1610//)i6(16 45//)i6, (17) 
We denote the vev of the SU{5) components of the 10//, containing all the possible SU{2) doublets, as 

05 = (5sc/(5)), 05 = (5sc/(5))- (18) 
If 05 = this operator contributes only in Y,, and Yd as usual in SU(5). In fact from Appendix A we have 

Mu - h'a5{bi-4b24) + h''^a5{bi+b2i), (19) 

M^ = 5h'a5bi + h'^a5{-3bi-3b24), (20) 

Md = h'a^{-3bi + 2b24) + h'^a^{bi+b24), (21) 

MJ = h'a^{^3bi~3b24) + h''^a^{bi+6b24). (22) 



2 This is the extra U{1) contained in SO{10) D SU{5) X UxW- 
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Again, because of the 45jy, the resuhing mass matrix is not symmetric but a generic matrix and it can contribute 
to the down-quark and lepton masses. Note that with a 45^f in B-L direction the resuhing mass matrix is 
antisymmetric in fact 

Mu = h'a5{b3~bi5) + h''^a5bi5, (23) 
AU = 5h'a5{b3 + 3bi5)-h''^a53bi5, (24) 
Md = h'a^{-3b3-bi5) + h'^a^bi5, (25) 
Mj = h'a^{-3bi5) + h''^a^{~b3 + 3bi5). (26) 

and putting 63 = the mass matrices are clearly antisymmetric. 

• Adding an 50(10) scalar singlet 1h we can consider also the dimension five operators 

1616126//!^, 161610h1h, 1616120^1^, (27) 

that behave in the same way as the renormalizable ones (see appendix A) introduced in eq. ([2]). Also these 
operators can be obtained integrating out a 16;,; — 16^ of heavy messenger fermions 

(16T26jf)i6(161^^)ye, (1610jf)i6(161ff)i6, (16120^)16(161^)16. (28) 

We see that the key ingredient to obtain type-I seesaw and TBM mixing is that the up-type SU{2) Higgs doublets 
in the 5io and 5i2o do not have vevs and so that they are not in linear combination of the light Higgs doublet. This 
can be a potentially problem since 5io takes a vev and it is mixed with the other light Higgs doublets. However the 
study of the complete scalar potential is beyond the scope of this paper and will be studied elsewhere. 



III. MODELS FOR TBM AND FIT OF FERMION MASSES AND CKM 

In Ref. [13] has been studied a model for TBM mixing with dominant type-H seesaw mechanism given in eq. ([2]) . In 
this section we present some possible modifications of the model given in eq. (HJ . In the models we will present below, 
TBM arises from both type-I and type-II seesaw mechanisms differently from Ref.0, where dominant type-II 
seesaw mechanism has been assumed for neutrino masses. We remark that the main problem with type-I seesaw is 
that the tree-level operator 16 16 10 gives equal contribution to the up-quark and Dirac neutrino mass matrix. But 
in order to fit quark masses and mixings with TBM neutrino mixing, the structure of the two mass matrices must 
be very different, namely the up quark mass matrix must be hierarchical while the Dirac neutrino mass matrix must 
be of TBM-type as in eq. ^ or the identity. So we need to disentangle the two sectors, leaving Dirac and Majorana 
neutrino masses of TBM-type and Dirac charged fermions masses hierachical and almost diagonal. ^From the previous 
section it is clear that one possibility is to replace the operator 16 16 10 of eq. ^ with the operator 16 16 45 120. 

In the following we will assume an underlying parity, like in (isj . making all the mass matrices hermitian and so 
reducing the number of free parameters. Another way to reduce the sometimes high number of parameters is to assume 
that the 45 get vev in the B-L direction. In this case the fermion mass matrices are symmetric or antisymmetric and 
not arbitrary. 

Examples of models with TBM neutrino mixing are listed below. The details of the fit are given in appendix B and 
C. We fit all charged fermion masses, the two neutrino mass square differences, leptons and quarks mixings, and the 
CKM phase for a total of 18 observables. For the operators 1616120/^45// and 16 16 10// 45// we always take zero 
vev for the component 55/7(5) of 120// and 10//, as described in the previous section (05 = 0). 

• case A:w^ / 16 16 T6// 16// + /i 16 16 45// 120// + /i' 16 16 45// 10// 
where the 45// takes vev in a general direction, that is 61 and 624 (see eq. (H))), are both different from zero. The 
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mass matrices are : 

M„ = h + f, (29) 
Md = n[h{2p- - 3) + h^i2p- + 2) + h'], (30) 

024 024 

M. = Mce[M2^ + 7) + .-(2^ + 2)] + /.'(-5^) + /^-(4^)}, OD 

M.., = i/, (32) 

= Tfl/, (33) 
Ml = TLf, (34) 

where h and /i' are generic matrices, r,, Ce and 6, are combinations of vevs (see appendix A). 
Results: 

= 0.0050, (35) 

with 26 parameters. 

We note that with the 45^ taking vev in B-L direction the number of parameters is considerably reduced but 
a good fit can not be performed. 

case B:w = f 16 16T26ij 1^ + /i 16 16 45^ 120^ + h' 16 1645^ 10^ 

where the 45^ takes vev in B-L direction. The mass matrices are: 

M„ = h^ + r^f, (36) 

Md = n{h^ + h^ + f), (37) 

Me = ri(ce/i'^- 3/1^-3/), (38) 

= -3r2/, (39) 

Mr = rnf, (40) 

Ml = TLf. (41) 



Results: 



with 16 parameters (2 d.o.f). 



X' = 5.6, (42) 



case C: «; = / 16 16 126ij Ijj + /i 16 16 45^ 120ij + h' 16 16 120'^ Ir 
where the 45ij takes vev in B-L direction. The mass matrices are: 



Results: 



with 18 parameters. 



Mu = /i^ + ra/i^ + r2/, (43) 

Md = n{h^ + h^ + f), (44) 

Me = n{cfh^ + c^h^-3f), (45) 

M^o = -3r2/, (46) 

Mr = TRf, (47) 

Ml = TLf. (48) 



= 0.0015, (49) 
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The last case reproduces basically the same structure of the renormalizable case (eq. ([2])) with type-II seesaw 
dominance studied for example in ref. [l^ . with just one more parameter cf . We note that the analysis performed in 
[T^ is based on a previous set of data (before the T2K and MINOS recent results). For this reason we show also an 
updated fit for that case, that can be used for comparison: 

X^ = 0.14, di.T= 461863 (50) 

with 17 parameters (1 d.o.f), where dpT is a parameter introduced in ^123. We note that the goodness of the fit is 
substantially unchanged compared with the old analysis, showing that in this class of models it is possible to obtain 
the desired (very small before T2K or more sizeable now) corrections to zero 9i3 from the charged lepton sector, 
taking into account an appreciable amount of finetuning. In fact the neutrino mass matrix is of TBM-type and it is 
diagonalized by TBM mixing matrix. The charged fermion mass matrices have hierarchical structure. Assuming all 
the parameters to be real, the charged lepton mass matrix is diagonalized by a rotation matrix O' characterized by 
three angles ^[3, ^^'^ ^23- The lepton mixing matrix Viep is given by the product Viep = O'^ ■ Vtbm so we have 

(Viep)i3 ^(SIS - C13S12), (51) 

{Vlep)i2 = ;^(C12C13 + S12C13 + S13), (52) 
(Vie;,)23 ^ ;y|(-Ci2C23 + S23(C13 + S12S13)), (53) 

where Sij = sin and Cy = cos dlj . We can have a large value for the reactor angle in agreement with the result of 
the T2K collaboration, and at the same time {Viep)j^2 ~ l/%/3 and (V;ep)23 ^ 1/^2 fine-tuning the mixing angles 6*'^. 

We observe that from type-I and type-II seesaw mechanisms we have for all the cases presented above 

= Mi - M,o -^M^o ^{rL-—)f = r,f, (54) 

where we have used the fact that / = Note that only a combination of the r^^j^ parameters enters in the neutrino 
sector. So counting the number of free parameters, and rji are equivalent to one free parameter instead of two. 



IV. RENORMALIZABLE THEORY 



The dimension five operators assumed in the previous section can be obtained from a renormalizable theory inte- 
grating out heavy messengers fields. In general the operator 16 16 0a 0b can be obtained from 

u; = 16 16x0a + 16X6x06 + Mxl6xT6x (55) 

where x — X is a couple of sets of fermion messengers'^ and it gives rise to the operator 16 16 4>a (t>b at a scale E <C M^. 

Moreover it is easy to take a symmetry forbidding the direct tree level operators 16 16 0a, for 0a = IOh, 120/f, \2&h- 
For example we can take a Z2 symmetry acting as 

(16,0a,b)^ (16,-0a,6), (X,X)^-(X,X). (56) 



Below we report explicit examples of renormalizable models from which the effective dimension five superpotentials 
assumed in the previous section can be obtained: 

• case A 

The matter and scalar fields content of a possible renormalizable model that can give the effective superpotential 
of the case A is given by: 



If only one messenger field is assumed the effective Yukawa mass matrix is rank one. 



7 





16 


IOh 16h 45h 120ff 


Ix 16x 16x 45x 


^2 


+ 







then the renormahzable superpotential is 
w = 16 Ix T6i/ + 16 45x T6ff + 1616x120^ + 16 T6x45i/ + 16 16xl0i/ +M 16xT6x + A/il^lx + Af45 45x 45x. (57) 

• case B 

The matter and scalar field content of the model is 





16 


Ih IOh 45h 120^ 126h 


16x1 16x1 16x2 16x2 




+ 






Z'2 


+ 


+ - - - + 


+ + - - 



and the superpotential is given by 
w = 16 16x1 T26h + 16 16x1 Iff + 16 16x2 10^ + 16 16x2 120^ + 16X6x245^ + A'h 16xiT6xi + ^2 16x216x2 (58) 
• case C 

The matter and scalar field content of the model is 





16 


Iff 120ff 45h 120'^ 126ff 


16x1 16x1 16x2 16x2 


Z2 


+ 






Z'2 


+ 


+ - - + + 


+ + - - 



and the superpotential is given by 
w = 16 16x1 T26h + 16T6xi 1^ + 16 16xi 120^ + 16 16x2 120^ + 16 16x245 h + Mi 16xiT6xi + ^^2 16x216x2 • (59) 

V. CONCLUSIONS 

Neutrino mixing data are in well agreement with maximal atmospheric angle, tri-maximal solar angle and may 
be with a non-zero and quite large (namely of order of the Cabibbo angle) reactor angle. TBM mixing gives zero 
reactor angle however it can be a reasonable starting point. In fact in GUT framework large deviation of the 1 — 3 
angle can arise from the charged sector. However a simple picture for TBM in 5*0(10) is still missing. In order to 
approach the problem recently has been studied models where light-neutrino mass matrix arises only from type-II 
seesaw mechanism. In this paper we studied the possibility that both type-I and type-II seesaw mechanisms yield 
TBM neutrino mixing in a 5*0(10) model. We have assumed that the superpotential contains only dimension five 
non-renormalizable operators. We studied three different possible scenarios for TBM neutrino mixing In each case 
proposed we make the fits of all the fermion masses and mixing angle. One case corresponds to the model studied 
already in Ref. [l3l for type-II seesaw dominance, while the other two are new. 

We found in both cases a good fit of all the data including the recent T2K result. In particular for the first model 
we found an excellent fit (x^ = 0.005) but with a high number (26) of parameters. We therefore can conclude that 
this case can be considered as a good starting point for a flavour theory that can reduce the number of the free 
parameters of the theory (for example introducing a flavour symmetry). Moreover in this case we did not need to 
introduce extra 50(10) singlets. For the second case we obtained = 5.6 but with only 16 free parameters and 2 
degree of freedom, making this case the most predictive. For third case we found a very good fit ^ = 0.002 with 18 
free parameters. This case also can be considered as a good starting point for a complete flavour theory. Even if we 
needed to introduce one 5*0(10) singlet we consider the last two cases as the most promising for the moment. 

For the three cases proposed, we give possible renormahzable realizations where we have introduced messenger 
fields and extra Abelian symmetries. 

We remark that in this paper we focused on the flavour secotr and we do not make a full analysis of the model. 
In particular we leave to a future analysis the study of the Higgs potential and related issues such as the breaking 
pattern of 5*0(10) to the SM, problems related to the doublet-triplet splitting (proton-decay) and the achieving of 
exact coupling unification considering the possible breaking steps and the related threshold corrections. 
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Appendix A 

Here we report for convenience of the reader the table VIII of ref.[6]. In general an SO{10) dimension five operator 
can be written as 1616 0a</'b where (j)a.b are scalar fields 0a, & = Ijy, 16//, 16//, 45//,... and so on. For simplicity we 
assume that 50(10) is broken through SU(5) and we describe the contribution of the dimension five operators to 
the fermion mass matrices in the SU (5) language. When one of the components of 4>a and </>(, take vev and bi re- 
spectively (where i is the SU{5) index of the component), one generates contributions to the quark and lepton masses. 



case 


50(10) operator 


mass matrices 


IV 


(16Afl6//)io(16A/16H)io 


— Kn-h-, -)- K'^n^h- 


V 


(16a/T6//)i(16a/T6//)i 


= Ka^bi + K^aibs 
Ml = KsOsbs 

A/fjD — /<*" /71 /)1 


VI 


(16a/16//)45(16a/16//)45 


M — Kn^bT + n^b^\ 
Ml = -^Ksasbs 
Mr = -5Ksaibi 


VII 


(16AflO//)i6(16A/45j/)Te 


Mu = Kasibi - 4624) + ^'^05(61 + &24) 
M, = 5KarM + K'^ar,i-3bi - 8624) 
Md = Ka-g{-3bi + 2624) + K'^a-g{bi + 624) 
Mj = Ka^(-3bi - 3&24) + K'^a^ibi + 6624) 


VIII 


(16Afl20//)i6(16Af45i/)T6 


Mu = Ka45{bi - 4624) - K'^a45{bi + 624) 
M^ = 5Kar,bi - K'^as{-3bi - 8624) 

Md = K{a^ + aj^){-3bi + 2624) - K'^{a-^+ ajs){bi + 624) 
Mj = K{a^ - 3aj^){-3bi - 8624) - K'^ {a^ - 3aj^){bi + 6624) 


IX 


(16i\/16//)l2o(16A/16H)l20 


Md = K{a-^bi + 2ai%) + K'^ {aib^ + 2a-^bi) 
Mj = K{a^bi + 2ai%) + K'^{aib^ + 2a-gbi) 



TABLE I: The contributions to the mass matrices from S'O(10)-invariant dim-5 operators, from table VIII of ref.0|. K is an 
arbitrary matrix. 



Below we report the contributions to the mass matrices from 50(10) invariant renormalizable Yukawa couplings. 
Different VEVs of the same 50(10) Higgs multiplet carry a subscript indicating the SU{5) component they belong 
to. 



case 


50(10) operator 


mass matrices 


I 




= = Y10V5 
Md = M, = Vio% 


II 


16Ail6Afl20// 


Mu = yi20t'45 

M^ = Yl20«5 

Md = yi20(% + l'45) 

MT = Yi2o{v^ ~ 


III 


16Afl6AfT26// 


Mu = Yj^vs 
M„ = -3Yj2gV5 

Md = ^126^^ 

Me = -SYj^v^ 
Ml = Vt26 «i5 
Mr = Vt26 ^'i 
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Appendix B 

In this section we show the fitting procedure used in our analysis. For charged fermions and CKM mixings the fit 
are performed on the set of data evolved at the GUT scale showed in Tab. |lTl The threshold effects are not considered, 
because they are model dependent and we try to make a general analysis valid for the various models. In these theories 
there are no constrains on the value of tan(3, so we use the high scale evolved data in the case of tanP = 10. 



Observables 


Input data 


mu[MeV] 


0.55 ± 0.25 


mc[MeV] 


210 ± 21 


mt [GeV] 




nidlMeV] 


1.24 ±0.41 


ms[MeV] 


21.7 ±5.2 


mt[GeV] 


^^•""-0.09 


me[MeV] 


0.3585 ± 0.0003 


m^[MeV] 


75.672 ± 0.058 


m,r[GeV] 


1.2922 ± 0.0013 


V^s 


0.2243 ± 0.0016 




0.0351 ±0.0013 


Vub 


0.0032 ± 0.0005 


J X lO"'^ 


2.2 ±0.6 



TABLE II: GUT scale data for charged fermions for tgl3 = 10 (ref.[23],[23]) 

For neutrino masses and PMNS mixings we use the results in Tab. IIIII These values are obtained with a global fit 
considering also the recent results from T2K and MINOS. In the models we considered we never obtain degenerate 
neutrino mass spectrum, so the effects of the evolution from the low energy scale to the GUT scale can be considered 
negligible to a good approximation for these observables. 



Observable 


Input data 


Amii X 10~^[eU^] 


' •'-'^-0.18 


Amii X 10"^ [eU^] 


^•'J"-0.16 


sill? On 


013+0 007 

U.U10_Q Qog 




n 312+0017 

U.01Z_Q Qj^g 




^^2+0 06 



TABLE III: Neutrino masses and mixing in normal hierarchy (ref. 01) 
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Appendix C 

Here we give some other details on the results of the numerical analysis. In particular for the three cases we analysed 
we give the best fit parameters and the values for the observables that we obtain. 

• case A: 



h = 



h' = 



/ill 

hi2 e-'^" 



hi2 e 



i23 e 



l22 

-iSh 



ni 



'12 ' 



''12 ' 



h' 



22 



"is h 



23 ' 



^ /2 /l /l 

/l /2 + /o /l — /o 

V /i /i - /o h + h ) 



his e'^^is 
/las 

7 / i<5fc' 



'13' 



23 



e 23 



33 



(60) 



(61) 



(62) 



(a) 



Observable 


Best fit value 


mu[MeV] 


0.550 


mc[MeV] 


210 


mt[GeV] 


81.9 


mi[MeV] 


1.24 




21.7 


mt[GeV] 


1.06 


me[MeV] 


0.3585 


m^[MeV] 


75.67 


mr[GeV] 


1.292 


Vus 


0.224 


Vcb 


0.0351 


Vub 


0.00320 


J X lO"'^ 


2.20 


Amii X lQ-''[eV'^] 


7.59 


Amia X 10'^ [eV^] 


2.50 




0.0132 




0.312 




0.516 




0.00500 



(b) 





Best fit value 


hiiVu[GeV] 


1.40 


hi2Vu[GeV] 


-2.45 


c 


-1.39 


hi,Vu [GeV] 


13.1 




0.232 


h22Vu[GeV] 


5.10 


h23Vu[GeV] 


15.0 


<5/l23 


1.81 


h33Vu[GeV] 


79.1 


h[^v4GeV] 


-6.42 


h[2Vu[GeV] 


-6.22 




0.901 


h[3VulGeV] 


2.62 


ShU 


-0.652 


h'22'Vu [GeV] 


4.55 


h'23Vu [GeV] 


-31.4 




-1.06 


h[y,v^[GeV] 


27.3 


hvu[GeV] 


-3.23 


hvu[GeV] 


-0.155 


f2Vu[GeV] 


0.987 


r\ 1 tan P 


-0.00418 


Ce 


-0.619 


61 


-1.12 




1.93 


V^jVu X 10~^ 


0.00947 



TABLE IV: Fit result for the case A (26 parameters) described in Sect. 3 



' hii hi2 his 

hi2 /l22 ^23 
\ his ^23 ^33 

(0 ai2 ai3 ^ 
0'12 (723 
CIS —0-23 J 

' h fl fl 

fl /2 + /o fl — fo 
\ fl fl~ fo f2 + fo 



(a) 



(b) 



Observable 


Best fit value 


m4MeV] 


0.465 


mc[MeV] 


210 


mt[GeV] 


81.5 


rna[Mc\-] 


2.95 


m, \MeV\ 


23.2 


m6[GeV] 


1.09 


m^{MeV\ 


0.3585 




75.67 


mr[GeV] 


1.292 


Vus 


0.224 




0.0352 


Vub 


0.00321 


J X 10"'^ 


2.15 


Amil X 10"' [el^^] 


7.59 


Amiz X 10"^ [eF^] 


2.50 


sin^9i3 


0.0118 




0.315 


sin^ 023 


0.516 




5.64 



Parameter 


Best fit value 


hiiVulGeV] 


0.785 


hi2Vu[GeV] 


0.346 


h^,v„lGeV] 


7.83 


h22ih, [Oc\-] 


3.09 


h23Vu[GeV] 


4.35 


h33Vu[GeV] 


82.2 


CTl2Vn[GeV] 


-0.298 


TisVu [GeV] 


0.229 


(T23Vn [GeV] 


2.29 


}oVu[GeV] 


-0.935 


fiVu[GeV] 


0.173 


f2Vu[GeV] 


0.0316 


r\ 1 tan 


0.0132 


Ce 


1.15 


r2 


2.50 


V^/Vu X 10~^ 


0.0249 



TABLE V: Fit result for the case B (16 parameters) described in Sect. 3 
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• case C: 



/ /ill hi3 

= hi2 /l22 h23 
\ his ^23 ^33 

(Tl2 (713 
/l^ = « I (712 (T23 
CTl3 — ^23 

/2 /l /l 
/ = I /l /2 + /O /l - /O 
/l /l - /o /2 + /o 



(66) 



(67) 



(68) 



(a) 



Observable 


Best fit value 


m4MeV] 


0.550 


mc[MeV] 


210 


mt[GeV] 


82.2 


rridlMeV] 


1.24 


m,[MeV] 


21.6 


mi[GeV] 


1.06 


me[MeV] 


0.3585 


mfj,[MeV] 


75.67 


mr[GeV] 


1.292 




0.224 




0.0351 




0.00320 


J X 10"'' 


2.20 




7.59 


Amia X 10"^ [eF^] 


2.50 




0.0131 


sin^ 012 


0.312 


sin^ 023 


0.520 




0.00149 



(b) 



Parameter 


Best fit value 


hiiVulGeV] 


0.584 


hi2Vu[GeV] 


-0.548 


hi3Vu[GeV] 


-5.49 


h22Vu[GeV] 


3.55 


h23Vu[GeV] 


3.99 


h33Vu[GeV] 


81.8 


(^i2Vu[GeV] 


-0.317 


ai3Vu[GeV] 


2.79 


o-23«u [GeV] 


-7.09 


foVu[GeV] 


-0.999 


fiVulGeV] 


-0.207 


f2Vu[GeV] 


0.0290 


ri / tan /3 


0.0129 


r2 


1.85 


rs 


1.37 


cf 


1.48 




1.17 




0.0287 



TABLE VI: Fit result for the case C (18 parameters) described in Sect. 3 
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